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ABSTRACT 

The Kepler satellite has discovered a number of transiting planets around close binary stars. These 
circumbinary systems have highly aligned planetary and binary orbits. In this paper, we explore how 
the mutual inclination between the planetary and binary orbits may reflect the physical conditions of 
the assembly of protoplanetary disks and the interaction between protostellar binaries and circumbi- 
nary disks. Given the turbulent nature of star-forming molecular clouds, it is possible that the gas 
falling onto the outer region of a circumbinary disk and the central protostellar binary have different 
axes of rotation. Thus, the newly assembled circumbinary disk can be misaligned with respect to the 
binary. However, the gravitational torque from the binary produces a warp and twist in the disk, 
and the back-reaction torque tends to align the disk and the binary orbital plane. We present a 
new, analytic calculation of this alignment torque, and show that the binary-disk inclination angle 
can be reduced appreciably after the binary accretes a few percent of its mass from the disk. Our 
calculation suggests that in the absence of other disturbances, circumbinary disks and planets around 
close (sub-AU) stellar binaries, for which mass accretion onto the proto-binary is very likely to have 
occurred, are expected to be highly aligned with the binary orbits, while disks and planets around 
wide binaries can be misaligned. 

Subject headings: accretion, accretion disks - hydrodynamics - planetary systems: protoplanetary 
disks - stars: binary 



1. INTRODUCTION 

The extremely precise photometry and nearly contin- 
uous observations provided by the Kepler satellite have 
led to the discovery of a number of transiting planetary 
systems around stellar binaries. At the time of this writ- 
ing, six such circumbinary systems are known, includ- 
ing Kepler-16 (with stellar binary period of 41 days and 
planet orbital period of 229 days; Doyle et al. 2011), 
Kepler-34 (28 d, 289 d), Kepler-35 (21 d, 131 d; Welsh 
et al. 2012), Kepler-38 (19 d, 106 d; Orosz et al. 2012a), 
Kepler-47 (stellar binary orbit 7.45 d, with two planets 
of periods 49.5 d and 303.2 d; Orosz et al. 2012b), and 
KIC 4862625 (20 d, 138 d; Schwamb et al. 2012, Rostov 
et al. 2012). The stars in these systems have masses of 
order of the mass of the sun or smaller, and the planets 
have radii ranging from 3.0i?e (Kepler-47b) to 0.76i?j 
(Kepler-34b). 

By virtue of their detection methods, all the Kepler 
circumbinary systems have highly aligned planetary and 
stellar orbits, with the mutual orbital inclinations con- 
strained between 6 ~ 0.2° (for Kepler-38b) and 0^2° 
(Kepler-34b and Kepler-35b). In Kepler-16, measure- 
ment of the Rossiter-McLaughlin effect further indicates 
that the spin of the primary is aligned with the orbital 
angular momentum of the binary (Winn et al. 2011). 
A natural question arises: do misaligned (O J> 5°) cir- 
cumbinary planetary systems exist? If so, under what 
conditions can they form? 

One might expect that circumbinary systems naturally 
form with highly aligned orbits, since the associated or- 
bital angular momenta originate from the protostellar 
cores. However, several lines of evidence suggest that 
misaligned configurations may be present in some sys- 



tems: 

(i) Solar-type main-sequence binaries with large sepa- 
rations (J> 40 AU) often have a rotation axis misaligned 
relative to the orbital angular momentum (Hale 1994). 
Misalignments are also observed in some short-period bi- 
naries, such as DI Hercules (with orbital period of 10 
days; Albrecht et al. 2009; see also Albrecht et al. 2011; 
Konopacky et al. 2012; Triaud et al. 2012). 

(ii) Some binary young stellar objects (YSOs) are 
observed to contain circumstellar disks that are mis- 
aligned with the binary orbital plane (e.g., Stapelfeldt 
et al. 1998). Also, several unresolved YSOs or pre- main 
sequence binaries have jets along different directions, 
again suggesting misaligned disks (e.g., Davis, Mundt & 
Eisloffel 1994; Roccatagliata et al. 2011). 

(iii) Imaging of circumbinary debris disks show that 
the disk plane and the binary orbital plane are aligned for 
some systems (such as a CrB, f3 Tri and HD 98800), and 
misaligned for others (such as 99 Herculis, with mutual 
inclination > 30°; see Kennedy et al. 2012a, b). Also, the 
pre-main sequence binary KH 15D is surrounded by a 
precessing circumbinary disk inclined with respect to the 
binary plane by 10°-20° (e.g., Winn et al. 2004; Chiang 
k Murray-Clay 2004; Capelo et al. 2012), and the FS 
Tauri circumbinary disk appears to be misaligned with 
the circumstellar disk (Hioki et al. 2011). 

While the aforementioned "misalignments" may have 
various origins (e.g., dynamical interactions in few body 
systems), in this paper we focus on the possible existence 
of warped, misaligned disks around proto-stellar binaries. 
We consider scenarios for the assembly of circumbinary 
disks in the context of binary star formation (Section 
2). These scenarios suggest that circumbinary disks may 
form with misaligned orientations with respect to the bi- 
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nary. We then study the mutual gravitational interaction 
between the misaligned disk and the binary (Section 3) 
and the long-term evolution of the binary-disk systems 
(Section 4). We discuss our results in Section 5 and con- 
clude in Section 6. 

2. FORMATION OF BINARIES AND CIRCUMBINARY 
DISKS: SCENARIOS 

Binary stars are thought to form by fragmentation in- 
side the collapsing cores/clumps of molecular clouds, ei- 
ther due to turbulent fluctuation in the core ("turbu- 
lent fragmentation"; e.g., Goodwin et al. 2007; Offher et 
al. 2010) or due to gravitational instability in the result- 
ing disk ("disk fragmentation"; e.g., Adams et al. 1989; 
Kratter et al. 2008). In the turbulent fragmentation sce- 
nario, the binaries form earlier and have initial separa- 
tions of order 1000 AU. Disk fragmentation also leads 
to binaries with large initial separations (~ 100 AU). In 
both cases, continued mass accretion and inward migra- 
tion, either due to binary-disk interactions (e.g., Arty- 
mowicz & Lubow 1996) or dynamical interactions in few- 
body systems, are needed in order to produce close (sub- 
AU) binaries. Planet formation can take place in the 
circumbinary disk during or after the binary orbital de- 
cay. 

In the simplest picture, the proto-binary and circumbi- 
nary disk rotate in the same direction. However, molec- 
ular clouds and their collapsing cores are turbulent (see 
McKee & Ostriker 2007; Klessen 2011). It is natural that 
the condensing and accreting cores contain gas which ro- 
tates around different directions. Even if the cores are 
not turbulent, tidal torques between neighboring cores in 
a crowded star formation region can change the rotation 
direction of the outer regions of the condensing/ accreting 
cores. Thus the gas that falls onto the central protostel- 
lar core and assembles onto the disk at different times 
may rotate in different directions. Such "chaotic" star 
formation has been seen in some numerical simulations 
(Bate et al. 2003). In this scenario, it is reasonable to 
expect a rapidly rotating central proto-stellar core which 
fragments into a binary, surrounded by a misaligned cir- 
cumbinary disk which form as a result of continued gas 
accretion. 

The mutual gravitational interaction between a proto- 
binary and the circumbinary disk leads to secular evo- 
lution of the relative inclination between the disk and 
the binary plane. In most cases, this interaction, com- 
bined with continued mass accretion, tend to reduce the 
misalignment. We will address these issues in the next 
two sections. Note that previous works have focused on 
warped circumstellar disks inclined relative to the bi- 
nary (e.g., Papaloizou & Terquem 1995; Bate et al. 2000; 
Lubow & Ogilvie 2000). The warped/twisted circumbi- 
nary disks studied below have qualitatively different be- 
haviours. 

3. WARPED CIRCUMBINARY DISKS 

3.1. Disk-Binary Interaction 

Consider a circumbinary disk surrounding a stellar bi- 
nary. The two stars have masses Mi and M2, and are 
assumed to be on a circular orbit with semi-major axis 
a. The circumbinary disk has surface density E(r), and 
extends from r[ n to r out (3> ri n ). The inner disk is trun- 
cated by the tidal torque from the binary, and typically 



ri n ~ 2a (Artymowicz & Lubow 1994; MacFadyen & 
Milosavljevic 2008). The orientation of the disk at ra- 
dius r (from the center of mass of the binary) is specified 
by the unit normal vector l(r). Averaging over the bi- 
nary orbital period and the disk azimuthal direction, the 
binary imposes a torque per unit area on the disk element 
at radius r given, to leading order in a/r, by 



3 GM t r]T,a 2 



(h-l)(l b xl), 



(1) 



where M t = Mi + M 2 is the total mass, r\ = M 1 M 2 /M? 

the symmetric mass ratio of the binary, and h is the 
unit vector along the orbital angular momentum of the 
binary. Under the influence of this torque, the angular 
momentum of an isolated disk element would precess at 
the frequency — Q p cos ft, where j3 is the angle between lb 
and Z, and 

3ry a 2 



4 r 2 



n(r), 



(2) 



with Q(r) ~ tt K = (GMt/r 3 ) 1 / 2 the disk rotation rate. 
Since fi p depends on r, the differential precession can 
lead to the warping (change with r of the angle between 

I and h) and twisting (change of I orthogonal to the I — lb 
plane) of the disk. 

3.2. Dynamical Warp Equations for Low-Viscosity disks 

Theoretical studies of warped disks (Papaloizou & 
Pringle 1983; Papaloizou & Lin 1995) have shown that 
there are two dynamical regimes for the linear propaga- 
tion of warps in an accretion disk. For high viscosity 
Keplerian disks with a J; S = H/r (where H is the half- 
thickness of the disk, and a is the Shakura-Sunyaev pa- 
rameter such that the viscosity is v = aH 2 £l), the warp 
satisfies a diffusion-type equation with diffusion coeffi- 
cient 1/2 = v/(2a 2 ). For low-viscosity disks (a <^ 5), on 
the other hand, the warp propagates as bending waves at 
about half the sound speed, c s /2. Protoplanetary disks 
with a ~ 10 _4 -10 -2 likely satisfy a <, S (e.g., Terquem 
2008; Bai & Stone 2011). For such disks, the warp evolu- 
tion equations governing long-wavelength bending waves 
in the linear regime (\dl/dhxr\ <C 1) are given by (Lubow 
& Ogilvie 2000; see also Lubow et al. 2002, Ogilvie 2006) 



£r 2 ft 



dl 1 <9G 



dt r dr 

dG _ m 2 - a 2 

dt ~ { 2il 



I x G 



YflG 



Y*H 2 n 2 z r 3 n dl 
4 dr 



(3) 
(4) 



where fl r and VL Z are the radial epicyclic frequency and 
the vertical oscillation frequency, G is the internal torque 
of the disk, and the surface density S(r) is taken to be 
the same as that of the unwarped disk. These equations 
are only valid for a <, S < 1, |£! 2 - Q 2 \ <, D, 2 S and 
|f2 2 — f2 2 | <S fl 2 S. For circumbinary disks considered 
here, the rotation rate differs from the Keplerian rate 
by an amount \fl 2 - n 2 K \/n 2 K = O (na 2 /r 2 ) + O (S 2 ), 
and similarly for fl r and fi z (see Eqs (|17 | - (fT8|l below). 
So the validity of equations ([3]) -((4]) requires rja 2 jr 2 <, 6, 
a condition that is generally satisfied. 
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Fig. 1. — Steady-state warp (top panel) and twist (bottom panel) 
profiles, in degrees, of a Keplerian circumbinary disk for which 
the outer disk is misaligned by 10° with respect to the angular 
momentum axis of the binary. Both stars in the binary have the 
same mass, and the disk parameters are p = 0.5, <5 = 0.1, r; n = 2a, 
and a = 0.01 or a = 0.001. The profiles for a = 0.001 are rescalcd 
to show the approximate scaling of the warp and twist with a. 

In the absence of the external torque (Tt = 0), equa- 
tions ©-(U]) admit wave solutions. If we define a Carte- 
sian coordinate system so that l z ~ 1 and \l x ,y\ *C 1, 
then a local (WKB) bending wave with l x ,y, G oc e lkr - luJt 
has a phase velocity cj/fc ~ ±c s /2 = ±H£l/2 (assuming 

LJ <C fl ~ fir — fla). 

3.3. Steady- State Warp and Twist of Circumbinary 

disks 

We now consider a circumbinary disk whose rotation 
axis at the outer radius (r out ), i ut = I (rout), is in- 
clined relative to the binary direction lb by a finite angle, 
/3(^out) = ©• This corresponds to the situation where the 
outer disk region is fed by gas rotating around the axis 
lout- In the Cartesian coordinate system with the z-axis 
along £&, the disk direction l(r) can be written as 

l(r) — (sin j3 cos 7, sin j3 sin 7, cos j3), (5) 

where /3(r) is the warp angle and "f(r) is the twist an- 
gle. At r = r out , we have Z ou t = (sin 6,0, cos <d) without 
loss of generality. A steady-state warp/twist is reached 
after a few bending wave propagation times across the 
whole disk. The steady-state warp/twist profile can be 
obtained by numerically integrating Eqs. ©-(U) and set- 
ting d/dt = 0. Figure 1 depicts selected numerical re- 
sults. 

Physically, the steady-state warp/twist profile is de- 
termined by balancing the internal viscous torque G of 
the disk and the external torque T^. The viscous damp- 
ing time of the disk warp [associated with the viscosity 
v 2 = vj(2a 2 )\ is 



tv2 — — 

V2 



2a 2 r 2 



2a 
SHI' 



(6) 



A critical warp radius r warp is obtained by equating t v2 
to the precession time fi" 1 of an isolated disk element, 
i.e. t v2 Q p = 1 at r = r 



where 

3an / a 
~2~ 



(7) 



This gives 



' warp 



( San\ 



1/2 



(8) 



For r warp 3> ri n , we expect that, in steady state, the disk 
well inside 



warp 



be aligned with the binary Z&, while 



the disk well outside 



warp 



be aligned with l c 



How- 
ever, if the inner disk radius r- ln is larger than r warp , or 
(t v2 Q p )- m 5s 1 (the subscript "in" means that the quan- 
tity is evaluated at r = r- m ), then the whole disk can be 

approximately aligned with Z out , with very small warp 
between the inner and outer edge of the disk. For stan- 
dard disk parameters (e.g., n ~ 1/4, a ~ 10~ 2 , S ~ 0.1, 
r in ~ 2a), the inequaility (t v2 £l p ) in <, 1 or 

^"warp 5s r[ n is 

well satisfied. 

For Keplerian disks (Q r = Q z = f2) and small warps, 
we can obtain approximate, analytic expressions for the 
disk warp and twist (see Foucart & Lai 2011 for a simi- 
lar calculation of magnetically driven warped disks) . For 
concreteness, we consider disk models with constant a, 
and assume that the surface density and the dimension- 
less disk thickness have the power-law profiles 



E(xr- p , 



-ocr^- 1 '/ 4 . 



(9) 



so that M ~ vY, = aH 2 QY =constant Q. Equations 
0-flU) then reduce to 



d_ 

Or 



3/2 



dr 



AarT h 



(<5 2 r 5 £fi 2 ) 



(10) 



We adopt the zero-torque boundary condition dl/dr — 
at the inner disk radius r = r, n . Since Tb falls off 
rapidly with r as T;, oc r~ 3-p , we can integrate Eq. (|10[) 
approximately to obtain: 



3/2 



d 



4a 



P 



5 2 r^V? 



1+p 



Using the outer boundary condition l(r out ) = lout, we 
find 



l(r)-i 



80 



T 6 (r) 



where 



l+p J (6 2 r 2 j:n 2 ) i: 



r \ p+1 1 
\r\nj 2p + 



3/2 

— ) F p (r) 



3' 



(12) 



(13) 



and t v2 Sl p is given by Eq. (J7J. In deriving Eq. (TT21) . we 
have assumed that the warp and twist of the disk are 
small compared to its inclination relative to the binary 
axis, i.e., \l(r) - l out \ <C sin/3, or (t v2 Q p ) in *C 1. The 

1 In practice, these scalings are unlikely to be valid for the entire 
radial extent of the disk. As the warp and twist of the disk are 
mostly due to the torque acting at small radii (Tj, oc r~ s ~ p ), we 
are only concerned about the approximate value of p for r close to 



4 



total change in I across the disk is then 

8(i v 2^p)in 



In 



2p + 3 



(lb ■ I out) (I b x lout)- 



The net twist angle across the disk, A7 tw i s t 
is 



7ir 



(14) 

■ 7out , 



A 7t , 



8(£ V 2^p 



cos 



(2p + 3) 
12 l arj 



(2p + 3) \5t 



cos 0. 



(15) 



For Keplerian disks, the warping is only a second-order 
effect: to first order ij n — Z out is perpendicular to Z ou t, and 
the disk is only twisted. The torque acting on the inner 
disk does however have a component in the lb— lout plane, 
due to the small difference in orientation between Zj n and 
l out . The net warp angle, A/3 warp = An - /3(r out ) = 
An - 0, is given by 



4(£ V 2^p)i 



2p 
6 



sin(20) 



2p + 3 



at] 

61 



sin(20). (16) 



As noted before, these expressions for A7t w ist and A/3 warp 
are valid only for |A7 twist | < 1, or (t v2 fl p )i n < 1. 

Figure 1 shows that the numerically integrated warp 
profile agree with both the analytic amplitude of the 
warp and its scaling with the viscosity parameter a. 
Thus, for standard disk and binary parameters (r/ = 0.25, 
a = 10~ 3 -10~ 2 , 5 — 0.1 and r; n ~ 2a), the steady-state 
disk is almost flat, with its orientation determined by 
i ou t, i.e., the angular momentum axis of the gas falling 
onto the outer disk. 

Deviations from a Keplerian disk profile can modify 
the above results, as differences between the epicyclic 
and orbital frequency of the disk are induced by both 
the finite thickness of the disk and the deviation of the 
binary gravitational potential from its point-mass value. 
To first order in 5 2 and r/(a/r) 2 , we have (assuming small 
binary-disk inclination 0) 



n 2 



GM t 



GM t 



1 + 



1 



3f]a 2 
4r 2 
3r]a 2 
4r 2 



CS 2 
D5 2 



(17) 
(18) 



where C and D are constants of order unity which 
depend on the density/pressure profile of the disk, 
and the epicyclic frequency was computed from fi 2 = 
(2n/r)d(r 2 n)/dr. We thus have 



n 2 -n 2 
n 2 



3-qa 2 
~2^ 



+ (D- C)S 2 



(19) 



It is worth noting that for 6 = constant, the S 2 term 
vanishes (since C = D in that case) . Including this result 
in Eq. (j4j leads to an additional warp 



A/3, 



NK 
warp 



-K±(±) sin(20) 



(20) 



with 



K 



0.97] 
2p+7 



ftS 2 



(21) 



and k a constant depending on the profile of 5(r) close 
to r; n (k — for constant S, and of order unity for slowly 
varying 5). Although this warp remains relatively small 
for most disk parameters (A/3^ p ~ 0.1/3), for a <S 0.03 
it will be larger than the second-order Keplerian warp 
given by Eq. (fl"6)) . As most of the torque on the disk 
is due to the contributions at radii r ~ rj n , this warp 
will also reduce the Keplerian change Al = I — l out [See 
Eq. ([T2jl ] in the orientation of the disk by a factor of 
order (sin An/ sin Aut)- 

4. EVOLUTION OF THE RELATIVE BINARY - DISK 
INCLINATION 

As discussed in Section 3, the binary torque [Eq. (JTJ] 
induces a small warp and twist in the circumbinary disk. 
The steady-state warp/twist is achieved on the timescale 
of bending wave propagation 



dr 



HQ/2 5 out n o 



(22) 



This timescale is much shorter than the age of the system 
or the gas accretion time 



1 



a^out^out 



(23) 



On a timescale longer than t waxp , the warped disk ex- 
erts a back-reaction torque on the binary, aligning lb 
with the disk axis (more precisely, with Z ou t). To de- 
termine this torque, recall that in the Cartesian coor- 
dinate system that we have set up, lb — (0, 0, 1) and 
'out = (sin 0, 0, cos0). In the small- warp approxima- 
tion, l(r) ~ (sin©, AL, cos 0), where Al y is the y- 
component of l(r) — l ut, given by [see Eq. (fT2]) ] 



AL 



4 t v 2^lp 



F p (r) cos sin 0. 



(24) 



Note that |AZ X | < \Al y \ < sin© when t v2 ft p < 1. Thus 
the torque on the disk element, Eq. ([!}, becomes 



3 GM t T]iZa 2 



■s© (- 



A; y ,sin©,0 



The back-reaction torque on the binary is 

/Tout 

T= - 2nrT b dr 

J Ha 

The a;-componnet of T tends to align lb with Z out : 



(25) 



(26) 



T x 



72-k 



(2p + 3)(4p + 5) V Si 



CX7] 



J a 4 E in ri 2 n cos 2 0sin0. 



( 27 ) 

If we define the angular momentum of the inner disk 
region by 



(AJ)i n EE 27T(I> 4 f2)i 



(28) 



5 



then T x can be written as 

r 32 

I X 



(2p + 3)(4p- 



— (AJ) in (^t v2 ) in , (29) 
5) 



where t v2 Cl p = (3ar]/2)(a/rS) 2 [see Eq. Q] 
The y-component of T on the binary is 



8 " GM tV ^^ cos6 sine. 



2(1 +p) 



(30) 



This makes the binary axis h precess around Z ou t at the 
rate 



3ir ( Sj n a 



2(1 +p) \M t r h 



fl b cosQl out , (31) 



where fib = (GMt/a 3 ) 1 / 2 is the orbital frequency of the 
binary. Since T y does not induce permenant change of 
the inclination angle Q, we will not consider it further in 
this paper. 

The expressions for T that we just derived are valid 
for Keplerian disks. It is worth noting, however, that the 
back-reaction torque is mostly independent of the non- 
Keplerian warp computed in Eq. (|20l) : to leading order, 

T x is proportional to the twist Al y and Ty to sine. The 
only effect of the non-Keplerian warp is thus to modify 
T by a factor of order (sin/3; n / sin/3 out ). 

Mass accretion from the circumbinary disk onto the 
binary can also contribute to the alignment torque. The 
accretion streams from r; n will likely land in both stars, 
probably through circumstellar disks (e.g., Artymowitz 
& Lubow 1994; MacFadyen & Milosavljevic 2008). Given 
the complexity of the process, we parametrize the align- 
ment torque due to accretion as 



% cc ,x =gM(GM t r in f/ 2 sinQ, 



(32) 



where g is a dimensionless number of order unity. In 
writing Eq. (|32|) . we have used the result of Section 3 
that the steady-state circumbinary disk is only slightly 
warped \J3(r- m ) ~ Q]. Since the mass accretion rate is 



given by M 
Eq. (H7J) as 



"SnvY, = 37ra((5 r f2£)j n , we can rewrite 



T x /M(GM t r in ) 1 / 2 cos 2 esine, 



where 



/ = 



24 



(33) 



(34) 



(2p + 3)(4p + 5) ' \8 in^ii 
The total alignment torque on the binary is then 

'/align — T&cCiX H~ Tx 

~ (g + / cos 2 Q)M(GM t r in ) 1/2 sin 6. (35) 

Assuming that the angular momentum of the binary, 
Lb = -qM^GMta) 1 / 2 , is much less than that of the disk 
(and the material falling onto the disk) , the torque 7Iii gn 
leads to alignment between lb and i ou t, on the timescale 



Malign 



Lb sin e f]M t ( a 



Tlign M V r m/ (,9 + /cOS 2 e) 



1/2 



(36) 



The secular evolution of 6(f) is determined by the equa- 
tion 

d e sine (37) 



dt 



Malign 



For 6 <C 1, this can be easily solved: Starting from the 
initial angle 0(ti), the inclination evolves according to 



e(t) = e(ii)exp 



AM /rinN 1 / 2 



r]M t \ a 



(t) '<•+/> 



(38) 



where AM is the total mass accreted through the disk 
during the time between t{ and t. 

5. DISCUSSION 

The calculations presented in Sections 3 and 4 show 
that a circumbinary disk formed with its rotation axis 
iout (at large distance) inclined with respect to the bi- 
nary angular momentum axis lb will attain a weakly 
warped/twisted state, such that the whole disk is nearly 

aligned with Z out (see Section 3). However, the inter- 
action torque between the disk and the binary tends to 
drive lb toward alignment with i out . The timescale of this 
alignment is given by Eq. (|36|) . and the relative binary- 
disk inclination 6 evolves according to Eq. (|38|) . 

Note that both the accretion torque and gravitational 
torque contribute to the alignment. If only the accre- 
tion torque were present (i.e., g ~ 1, / = 0), the align- 
ment timescale would be of the same order as the mass- 
doubling time of the binary (taiign ~ 4 x 10 7 yr for 
Mi = M 2 = 1 M , r in ~ 2a and M - 1O" 8 M yr" 1 ), 
and a significant fraction of the binary mass would have 
to be accreted (AM ~ 0.4Mq) in order to achieve an e- 
fold reduction of Q [see Eq. (pigj) ]. However, the gravita- 
tional torque dominates over the accretion torque, since 
the condition / 3> 1 can be satisfied for a wide range 
of disk/binary parameters (although a 2 f <C 1 must be 
satisifcd for our equations to be valid). For example, 
for p — 3/2 (the density index for the minimum solar 
nebula) , and r\ = 1/4 (equal mass binary) , we have 



(39) 



Thus, the alignment timescale is (for / 3> 1 and cos 2 Q 
1) 



^align — 2.5 



M 



yMt \ 

0.5MJ I lO- 8 M /yr, 



(40) 

(r)Mt is the reduced mass of the binary) . The amount of 
mass accretion needed for an e-fold reduction of Q is [see 
Eq. O] 



(41) 



(AM) e ^0.05(r?M t )(^) ^ 



Thus, only a small fraction of the binary mass has to be 
accreted to achieve significant reduction of 6. 

We comment on two assumptions adopted in our cal- 
culations of Sections 3-4: 

(i) We have assumed that the binary separation a is 
constant. In reality, the binary-disk alignment can take 
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place simultaneously as the binary orbit shrinks (due to 
binary-disk interactions). However, since the alignment 
timescale depends only on the ratio r- ln /a, we expect our 
result to be largely unchanged in such a situation as long 
as n n keeps track of a during the orbit decays. 

(ii) We have assumed that there is a constant supply of 
gas at the outer disk and the total angular momentum of 
the disk, Ldisk, is much larger than that of the binary, Lb- 
If we consider an isolated circumbinary disk (e.g., when 
an episode of mass infall from the turbulent cloud/core 
onto the central binary occurs) with Ldisk comparable 
or smaller than Lb, then both the binary and the disk 
will precess around a common total angular momentum 
axis while the mutual inclination © evolves. In this case, 
equation (l37t is replaced by 

dQ %i ign %n gn sin6 / L fc cos6\ 
— = - — — cos 8 = 1 H . 

dt Lb -^disk ^align \ -^disk / 

(42) 

This equation neglects the accretion torque, so that 
7Iii gn — l~x- Note that in general, T x and i a iign will be 
modified from the expressions given in Section 3. But 
as long as r out 3> r- m , we expect the corrections to be 
small (Foucart & Lai, in preparation). It is also of inter- 
est to note that when cos0 < — L^k/Lb, the gravita- 
tional torque tends to drive O toward 180° (i.e., counter- 
alignment). 

6. CONCLUSIONS AND IMPLICATIONS 

In this paper, we have considered scenarios for the 
assembly of proto-planetary disks around newly formed 
stellar binaries. The shape and inclination of the disk rel- 
ative to the binary will determine the orbital orientations 
of the circumbinary planets that are formed in the disk. 
Because of the turbulence in molecular clouds and dense 
cores, inside which protostellar binaries and circumbi- 
nary disks form, and also because of the tidal torques 
between nearby cores, it is possible, and even likely, that 
gas falling onto the outer region of the circumbinary disk 
rotates along a direction different from the rotation axis 
of the binary. Thus in general, the newly assembled cir- 
cumbinary disk will be misaligned with respect to the bi- 
nary. However, the gravitational torque from the binary 
produces a warp and twist in the disk, and the back- 
reaction torque associated with that twist tends to align 
(under most conditions) the disk and the binary orbital 
plane. We have presented new calculations of the inter- 
action between the warped/twisted disk and the binary, 
and showed that the disk warp is small under typical con- 
ditions. More importantly, we have derived new analytic 
expressions for the binary-disk alignment torque and the 
associated timescale [see Eq. (j40|)]. Our results show that 



the misalignment angle can be reduced appreciably after 
the binary accretes a few precent of its reduced mass [see 
Eq.|4l]. 

Proto-binaries formed by fragmentation (either turbu- 
lent fragmentation or disk fragmentation; see Section 2) 
have initial separations much larger than 1 AU. Signifi- 
cant inward migration must occur to produce close (sub- 
AU) binaries. Since mass accretion necessarily takes 
place during disk-driven binary migration, our results 
then suggest that close binaries are likely to have aligned 
circumbinary disks, while wider binaries can have mis- 
aligned disks. This can be tested by future observations. 
The circumbinary planetary systems discovered by Ke- 
pler (see Section 1) all contain close (period <, 41 d) 
binaries. If the planets form in the late phase of the 
circumbinary disk (as is likely to be the case considering 
the relatively small planet masses in the Kepler systems), 
then the planetary orbits will be highly aligned with the 
binary orbits, even if the initial disk has a misaligned 
orientation. This is indeed what is observed. 

Of course, given the complexity of the various processes 
involved, one may see some exceptions. In particular, in 
this paper we have not considered any dynamical pro- 
cesses (few body interactions) that may take place after 
the binary and planet formation. Such processes can also 
affect the mutual inclinations of circumbinary planets. 

Observationally, tertiary bodies on misaligned orbits 
around close, eclipsing binaries can be detected by 
searching for periodic eclipse timing variations. This has 
led to the identification of many binaries with tertiary 
companions (e.g., Liao & Qian 2010; Gies et al. 2012). 
A number of post-main-sequence eclipsing binaries have 
been claimed to host candidate circumbinary planets, 
such as HW Virginis (Lee et al. 2009), HU Aquarii 
(Qian et al. 2011; Gozdziewski et al. 2012), NN Ser- 
pentis (Beuermann et al. 2010), DP Leonis (Beuermann 
et al. 2011) and NY Vir (Qian et al. 2012). However, 
some of these claims are controversial since the proposed 
planetary orbits may be dynamically unstable on short 
timescales (see Horner et al. 2011,2012). Currently, no 
misaligned (inclination ^5°) circumbinary planets have 
been confirmed around main-sequence binaries. 

Overall, our calculations in this paper illustrate that 
the mutual inclinations and other orbital characteristics 
of circumbinary planetary systems can serve as a diagnos- 
tic tool for the assembly and evolution of protoplanetary 
disks and the condition of formation of these planetary 
systems. 
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